Circular geodesies and accretion disks in Janis-Newman-Winicour and Gamma metric 
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We study here circular timelike geodesies in the Janis-Newman-Winicour and Gamma metric 
spacetimes which contain a strong curvature naked singularity and reduce to the Schwarzschild 
metric for a specific value of one of the parameters. We show that for both the metrics the range 
of allowed parameters can be divided into three regimes where structure of the circular geodesies 
is qualitatively different. It follows that the properties of the accretion disks around such naked 
singularities can be significantly different from those of disks around black holes. This adds to 
previous studies showing that if naked singularities exist in nature, their observational signature 
would be significantly different from that of the black hole. 
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I. INTRODUCTION 

One of the most important unresolved problems in gen- 
eral relativity, and an issue that has far reaching obser- 
vational implications, is that of the final fate of complete 
gravitational collapse of a massive body such as a star. 

The Israel-Carter conjecture asserts that a Kerr black 
hole, which is described by two parameters, namely mass 
and angular momentum, would be formed as the endstate 
of any generic complete gravitational collapse and thus 
all other information regarding the nature of the mat- 
ter fields, symmetries, initial conditions will be radiated 
away in the process of collapse. The above statement re- 
mains at the stage of conjecture as of now because prov- 
ing it analytically or numerically is turning out to be a 
very difficult task. We also do not understand very well 
how matter would behave at energies beyond nuclear en- 
ergy scales all the way up to Planck scale as we have no 
theoretical understanding of evolution of the equation of 
state of matter at such energies. Thus the theoretical 
studies do not provide a definite answer to the question 
whether or not black holes are the only possible outcomes 
of the gravitational collapse. 

In fact in the last few years new theoretical models 
and observations seem to suggest that the true nature of 
the process could be more complicated. Investigations 
of gravitational collapse of various matter configurations 
in the context of the general theory of relativity over 
last couple of decades has revealed that the endstates 
of complete gravitational collapse could either be black 
holes or naked singularities [l|. 

Therefore it is natural to ask the question of what hap- 
pens during a realistic collapse of a massive body and 
what will be the endstate of such a collapse. The possi- 
bilities can be divided into three main categories: 
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- The body radiates away all higher multipole mo- 
ments and all asymmetry, thus forming a Kerr 
black hole. 

- Collapse will halt before all matter is squeezed into 
a spacetime singularity, thus creating a finite sized 
final object that is different from a black hole. 

- Matter will fall into the final spacetime singularity 
preserving its symmetry structure, thus creating a 
final configuration in the form of vacuum spacetime 
with a singularity (that can be covered or naked). 

The analytical solutions describing dynamical collapse 
of a rotating body away from spherical symmetry are few 
and not well understood at present. Therefore a possi- 
ble mechanism to account for the first hypothesis is still 
missing. For this reason it seems natural to turn atten- 
tion to the other two options. Ideas following the sec- 
ond hypothesis have been proposed |2( , but the nature of 
these objects and the physical processes that could lead 
to the formation of the same are not well understood. 
Furthermore recently some new scenarios have been pro- 
posed that include the possibility that gravitational col- 
lapse can asymptotically halt, leading to the formation of 
static configuration of matter that may or may not con- 
tain a naked singularity and which could be either finite 
or infinite in extent |3J . 

On the other hand, if black holes and naked singular- 
ities which are hypothetical astrophysical objects occur 
in nature, a question would be how they would be ob- 
servationally different. Therefore some researchers have 
started looking into the observational features that would 
distinguish these entities. The black holes and naked sin- 
gularities could have rather different properties and this 
could also possibly shed a light on the nature of the ex- 
isting sources, like the supermassive dark objects that 
dwell at the center of galaxies. Recent studies of gravi- 
tational lensing (see for example [4|) and accretion disks 
(see 0)0) Si an d 01) have brought out interesting charac- 
teristic differences for these objects which could possibly 
distinguish them from each other. Recently it was also 
suggested by some of us that ultrahigh energy collisions 
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and fluxes of the escaping collision products could be 
used for this purpose in certain situations |8j. 

In the present paper we investigate the properties of 
circular geodesies and accretion disks in two spacetimes 
that arise naturally while generalizing the Schwarzschild 
solution. One is the Janis-Newman-Winicour (JNW) so- 
lution and the other one is the so called Gamma metric 
(7- metric henceforth) . 

The Schwarzschild metric is the unique vacuum so- 
lution to Einstein equations under spherical symmetry 
which is asymptotically flat. Different metrics can be ob- 
tained from the Schwarzschild solution by relaxing one 
or more of the assumptions above. Firstly, relaxing the 
spherical symmetry, we can consider either a static defor- 
mation keeping the axial symmetry intact, or we could 
invoke a spin. One of the solutions that can be obtained 
in the former case is the 7- metric @ , whereas the space- 
time in the latter case is described by the Kerr metric. 
On the other hand the assumption of the empty space- 
time can be relaxed by invoking the presence of either an 
electromagnetic field or a scalar field. In the former case 
one obtains the Reissner-Nordstrom solution, whereas in 
the latter case, by invoking a massless scalar field, the 
JNW solution is obtained |10| . 

The uniqueness theorems give a privileged status to the 
Kerr and Schwarzschild spacetimes. Nevertheless since 
we do not know if in a realisticphysical scenario those 
will necessary be the resulting metrics it is of interest 
to study spacetimes such as JNW and the 7-metric that 
depart from the black hole models and see if they can 
be observationally different from their black hole coun- 
terparts. The main point that makes the 7-metric and 
the JNW-metric more appealing is that they are contin- 
uously connected to the Schwarzschild solution via the 
value of one parameter. By reducing this parameter to 
a specific value we recover exactly the Schwarzschild so- 
lution. The interesting point is that the Schwarzschild 
metric has a black hole whereas these two generaliza- 
tions contain naked singularities and therefore we can 
directly compare the results obtained for models of ac- 
cretion disks in these metrics with the Schwarzschildcase. 
Furthermore the parameters appearing in the JNW and 
7-metric have a clear and direct physical interpretation, 
as one represents the addition of a scalar field while the 
other represents an axial deformation. 

The first step towards the study of accretion disks is 
an understanding of the circular geodesies in the equa- 
torial plane of the spacetime. The circular geodesies 
in Kerr and Reissner-Nordstrom spacetimes containing 
black hole as well as the naked singularities depending on 
the values of the spin and charge parameters respectively 
was recently carried out in Q. In this paper we investi- 
gate the structure of circular geodesies for the JNW and 
7 spacetimes. 

We find that the properties of accretion disks can be 
divided in three regimes according to the values assumed 
by a parameter in the solutions. In one regime, the stable 
circular geodesies extend all the way up to singularity, in 



the second the stable circular orbits can exist in two dis- 
joint disks separated by a gap, one disk extending from 
singularity to a finite radius and another extending from 
a larger radius to infinity. Finally, in the third case a 
photon sphere surrounds the singularity, and stable cir- 
cular geodesies can exist from a particular radius above 
photon sphere to infinity. 



II. CIRCULAR EQUATORIAL GEODESICS IN 
THE JNW METRIC 

In this section we investigate the structure of the cir- 
cular geodesies in the JNW spacetime. As mentioned 
earlier the JNW metric is obtained as an extension of 
the Schwarzschild spacetime when a massless scalar field 
is present. In a naive sense, the black hole event hori- 
zon, i.e. a coordinate singularity in the Schwarzschild 
spacetime transforms to a real physical singularity which 
is naked. The JNW metric depends upon two parame- 
ters, one related to the mass and an additional parameter 
which can be interpreted as the strength of the scalar field 
or 'scalar charge'. 

The JNW metric can be written in the following way 



ds 2 = -Adt 2 + A- X dr 2 + Bdn 2 



(1) 



where dil 2 = dfl 2 + sin 2 6d<h 2 is the line element of a unit 



two-sphere and the functions A{r) and B(r) are given by 
the following expressions 



A(r) 
B(r) 



2r-r (/i-l) 
2r + r (^ + l) 

l [2r + r (/i + l)]* +1 



(2) 
(3) 



The two parameters mentioned above are r (that is re- 
lated to the mass) and the scalar charge fi. The range of 
\x is given by 11 € (1, 00) and he scalar field is given by 
the following expression, 



a 1 

ip = — m 

A* 



2r-r Q (p-l) 



2r + r (> + l) 



(4) 



where /i and a are linked by the relation /1 = 
^1 + Stt^t^J . The JNW spacetime contains a strong cur- 
vature singularity which is globally naked at 



1 , 

r s ing = -r (/U - 1) , 



(5) 



for /i 7^ 1. 

The Schwarzschild solution in the usual form can be 
recovered from the JNW solution when /i = 1 after a 
simple coordinate transformation of the radial coordinate 
given by f = r + r , with r = 2M which is the usual 
Schwarzschild radius of the event horizon (note that the 
radial coordinate r of the JNW metric when = 1 does 
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not correspond to the usual Schwarzschild coordinate f), 
and where M is related to the total mass of the central 
body as measured by observers at infinity. Therefore the 
case where /i is close to 1 corresponds to the small devia- 
tions from the Schwarzschild metric due to the presence 
of a small massless scalar field, whereas large values of \x 
correspond to the large deviations from the Schwarzschild 
geometry 



11] 



A. Geodesies in JNW spacetime: Basic equations 

We now investigate the structure of timelike equato- 
rial geodesies in the JNW spacetime. The JNW metric 
is spherically symmetric and static and therefore admits 
the two Killing vectors £ t = d t and = d$. As a conse- 
quence we have the following conserved quantities along 
geodesic motion. 



E 



A(r)U l 



L = g ap ^u p = B{r)U* , 



(6) 
(7) 



where U a = (U\ U r , U\ V*) is the four-velocity of the 
test particle. The constants of motion E and L can be 
interpreted as the conserved energy and orbital angular 
momentum per unit mass of the particle respectively. 

Spherical symmetry implies that the motion of the par- 
ticle will be restricted to a plane which can be chosen by 
gauge freedom to be the equatorial plane. Thus we set 
6 = | and U B = 0. 

From equations ((5]) and (J7J using the normalization 
condition for the four- velocity U a U a = — 1, we obtain 
the equation for the motion of a particle in the radial 
direction. 



f 2 



V 2 = E 2 



(8) 



Here V can be thought of as an effective potential for the 
test particle and it is given by a following expression, 



V = A(r) 1 



L 2 
B(r) 



(9) 



B. Circular geodesies 

From the analysis of the above effective potential we 
can obtain the structure of circular geodesies. A particle 
orbiting along a circular geodesies must obey the condi- 
tions 



U r = U r = , 



(10) 



as r must stay constant. Here the dot denotes a derivative 
with respect to the affine parameter for geodesic i.e. it 
is the proper time. 

From equations jSJ and (fTO]) it follows that for circular 
geodesies 



V = E , and 



dV 
dr 







(11) 



Further one needs to ensure that the circular orbit is sta- 
ble against the small perturbations in the radial direction 
and this implies that for a stable orbit the effective po- 
tential must admit a minimum: 



d 2 V 
dr 2 



> . 



(12) 



When the effective potential admits a maximum i.e. 
when < the orbit is unstable against small per- 
turbations. The condition for the marginally stable orbit 



is 



d 2 V 
dr 2 



0, corresponding to the inflection point for the 
effective potential. 

Imposing the conditions (fTT|) on (|9|) we get the ex- 
pression for the conserved energy and momentum of the 
particle orbiting along the geodesic of radius r as 



L 2 = 



ro [2r + r (^ + l)]^ +1 
4(2r-r ) [ 2r - r (p - l)]* -1 ' 



E 2 = 



2r 



2r — ro 



2r-r o 0u-l) 



2r + r (/i+ 1) 



from which it follows that we must have 



(13) 
(14) 

(15) 



for the energy and angular momentum to be real. Thus 
the circular geodesies can exists only for those values of 
r for which this condition holds. The energy and angular 
momentum diverge at the radius 



r p h 



ro 
2 



(16) 



which represents the radius at which photons move on cir- 
cular trajectories and it is known as photon sphere. The 
photon sphere plays a crucial role in the study of grav- 
itational lensing since its presence delimits a boundary 
below which lensing effects cannot be observed, there- 
fore implying that two spacetimes that possess a photon 
sphere at the same radius cannot be distinguished by the 
observation of lensing effects even though their compact 
sources might be different. 

Imposing the condition for the stability of the circular 
geodesies against small perturbations along the radial di- 
rection in equation ([9]) we get 



4r 2 - 8rr a + r 2 (/i 2 - 1) >0 



(17) 



which implies that the radius of the stable orbit must 
satisfy the following condition 



where 



r > r+ or r < r_ 



r± =r (l±^V5-V 2 ) 



(18) 



(19) 



The marginally stable orbits are located at these two val- 
ues of the radii. 
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FIG. 1: The value of the radial coordinate (in units of ^xf) 
for the singularity (thick line), photon sphere (thin line), r+ 
(dotted and dashed line) and r_ (dashed line) as functions of 
the parameter [i. It is clear that the structure of the circular 
geodesies will be qualitatively different in the three different 
ranges of the parameter, namely (l,2),(2,v5) and (V5,oo). 

In Fig. [1] are shown the values of the radial coordi- 
nates for the singularity, photon sphere, r + and r_ de- 
pending upon the parameter /x. Three distinct ranges 
of the parameters emerge where the structure of the 
circular geodesies will be qualitatively different, namely 
(i G (1,2), ju G (2, VS), (J, € (Vb,oo). We discuss these 
cases separately. 

The conserved energy and angular momentum for test 
particles as a function of the radius of the circular orbit 
are plotted in Fig. [2jind Fig[3] respectively, for different 
values of the parameter /i. The qualitative behaviour of 
these functions is also manifestly different in the three 
regimes of the parameter mentioned above. 

1. Case fie (1,2) 

For this range of values of parameters we have a fol- 
lowing relation, as seen from Fig. [T] 

r+ > r ph > r smg > r_ . (20) 

This implies that r_ does not exist and a photon sphere 
covers the singularity. The stable circular orbits can exist 
in the region given by, 

re(r + ,oo), (21) 

which is above the photon sphere. Unstable circular or- 
bits can exist between the photon sphere and r + and no 
circular orbits can exists below the photon sphere up to 
the singularity. 

In Fig|4] the general behaviour of the effective poten- 
tial is shown as a function of the radius for a fixed value 
of the parameter fi in the range considered and for dif- 
ferent values of the angular momentum L. The effective 
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FIG. 2: The conserved energy per unit mass E is plotted 
against the radius of the circular orbit r (in units of 537) 
for different values of \i. The curves located higher in the 
diagram correspond to lower values of fi and vice versa. For 
values of /j, G (1, 2), E decreases with the decreasing r, attains 
a minimum and again increases to become infinite before the 
singularity is reached. For values of fj, G (2, \/5), E initially 
decreases with decreasing r, attains a minimum followed by 
a maximum and decreases to a zero value at the singularity. 
Whereas for large values of /j, £ (yo, 00), E is monotonically 
decreasing with decreasing r and takes a zero value at the 
singularity. 

potential admits one maximum and one minimum which 
correspond to the unstable and stable circular orbits, re- 
spectively. The unstable circular orbits lie in the region 
between the photon sphere and r + . 

Note that both E(r) and L(r) become infinite at a cer- 
tain finite radius r in this case (as can be seen from Figs. 
[5] and . As we approach the singularity from infinity, 
the stable circular geodesies can exist up to the mini- 
mum of energy and angular momentum. The unstable 
orbits exist in the region between the minimum and the 
radial value where E(r) and L{r) blow up, which is the 
photon sphere. No circular orbits can exist below the 
photon sphere a result that is consistent with the earlier 
discussion. 

2. Case n G (2, y/Z) 

In this range of parameters we get the following rela- 
tion 

r + > r_ > r sing > r ph . (22) 

This implies that the photon sphere does not exist and 
there are two distinct regions where circular geodesies 
can exist namely, 

r e (r+,00) and r G (r smg ,r_) . (23) 

A new interesting feature arises in this case since the 
second region extends all the way up to singularity. Par- 
ticles reaching the innermost boundary of the exterior 
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FIG. 3: The conserved angular momentum per unit mass L 
is plotted against the radius of the circular orbit r (in units 
of 27jj ) for different values of fi. The curves located higher in 
the diagram correspond to lower values of n and vice versa. 
For values of /x £ (1,2), L has a minimum at a finite r and 
increases diverging as the singularity is approached and as r 
goes to infinity. For values of /x £ (2, V5) L vanishes at the 
singularity and has a local maximum and a local minimum 
before increasing indefinitely as r goes to infinity. Whereas 
for large values of fi £ (\/5, oo), L increases monotonically 
from zero at the singularity to infinity as r goes to infinity. 

V 




FIG. 4: The effective potential for radial motion V is plotted 
as a function of radial coordinate for fi = 1.5 £ (1,2) for 
different values of L. It can be seen that the V admits one 
minimum and one maximum, which corresponds to stable and 
unstable circular orbits respectively. 



accretion disk would plunge in the region of unstable or- 
bits to reach r_ where they would shock and circularize 
again, inspiralling onto the singularity. 

In FigJS] the general behaviour of the effective poten- 
tial for this range of the parameter is shown for different 
values of L. The effective potential admits two minima 
which correspond to the two stable circular orbits in the 
two regions we discussed above. There is also one lo- 
cal maximum between the two minima, which lies in the 
unstable region between the two stable ones. 

The qualitative behaviour of the functions E(r) and 
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FIG. 5: The effective potential for radial motion V is plotted 
as a function of radial coordinate =?j for /1 = 2.1 £ (2, %/5) 
for different values of L. It can be seen that V admits two 
minima and one maximum in between, which indicate the 
presence of two regions with stable circular orbits separated 
by one region of unstable circular orbits. 

L(r) in this range of parameter shows that both quan- 
tities do not blow up, instead as r decreases they reach 
a maximum and then decrease to zero as singularity is 
approached (see Figs. [2] and [3]). The existence of two 
regions where stable circular geodesies can exist can be 
derived also from the analysis of E(r) and L(r) as the 
two regions correspond to the portions where the two 
functions decrease with the decreasing r, namely from 
infinity to the first minimum and from the maximum to 
the singularity. Since the two functions do not blow up 
at any allowed radius, no photon sphere is present in the 
spacetimc. 

3. Case n £ (v5 ; 00 ) 

Within this range of values for the parameter r + and 
r_ do not exist as the solutions of equation (fTTJ|) are not 
real. For this range of parameters we have 

> r ph , (24) 

from which we see that the photon sphere is not present 
as well. This implies that the stable circular orbits can 
exist for 

r £ (r smg ,oo) , (25) 

i.e. from the singularity to infinity everywhere on the 
equatorial plane. 

In Fig. 6 the general behaviour of the effective potential 
for this range of parameter values for various angular mo- 
menta is shown. The effective potential admits only one 
minimum which indicates the presence of stable circular 
orbit. Since the potential does not admit maximum, un- 
stable circular orbits do not exists. The circular orbits 
exist all the way up to the singularity. 
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FIG. 6: The effective potential for radial motion V is plotted 
against r (in units of ^7) for fi — 3 6 (V5, 00) for different 
values of L. The effective potential admits only one mini- 
mum. Unstable orbits do not exists since V does not admit 
maximum. 



In this case the qualitative behaviour of the functions 
E(r) and L(r) is decreasing from infinity until they be- 
come zero at the singularity (see Figs. [2]and|3|). This 
implies that the stable circular orbits would exist all the 
way up to singularity, which is consistent with the above 
discussion. 



by the following expressions 
2M X 7 



F=\\- 

G- 
H 



2Mr 



r 2 - 2Mr + M 2 sin 2 9 
(r 2 - 2M r)^ 2 
[r 2 ~ 2Mr - M 2 sin 2 9)i 2 ~ 1 



(27) 
(28) 
(29) 



The total ADM mass as measured by an observer at in- 
finity associated with the 7-spacetime is given by M tot = 
7M and the range of the deformation parameter is over 
all the positive real numbers with the Schwarzschild 
spacetime recovered when 7 = 1, whereas 7 > 1 (7 < 1) 
corresponds to an oblate (prolate) geometry respectively. 
In the limit 7 = a flat spacetime is obtained. 

There is a strong curvature singularity in the 7 space- 
time whenever 7^1 and it is located at the radial value 



2M 



(30) 



The singularity is visible to all observers if 7 < 1 while 
for 7 > 1 it has a directional behaviour being visible only 
to observers not on the polar plane. So if we confine our- 
selves to the equatorial plane, the singularity is a naked 
singularity [T2I ] . 



III. CIRCULAR GEODESICS IN THE GAMMA 
METRIC 



A. Geodesies in 7-spacetime: basic equations 



In this section we study the structure of the circular 
geodesies in the 7 spacetime. The 7-metric is obtained 
by extending the Schwarzschild spacetime when spheri- 
cal symmetry is deformed to develop a prolate or oblate 
spheroid. 

The 7-metric is a two parameter family of spacetimes 
belonging to the Weyl class of static, axially symmetric, 
vacuum solutions of Einstein equations that are asymp- 
totically flat. In principle all the metrics belonging to 
the Weyl class are known since there is a one to one cor- 
respondence between such spacetimes and solutions of 
the Laplace equation in flat two-dimensional space. The 
interesting feature of the 7-metric is that it is continu- 
ously linked to the Schwarzschild metric via one of the 
parameters. The two parameters can be defined as a 
mass parameter M and the parameter 7, that quantifies 
the deformation from spherical symmetry. We investi- 
gate the circular geodesies in 7-spacetime as a function 
of the deformation parameter. 

The 7-metric is given by 



ds 2 



-Fdt 2 + F^lGdr 2 + Hd9 2 + (r 2 - 2mr) sin 5 



(26) 



where the functions F(r), G(r,0) and H(r,9) are given 



The complete structure of geodesic motion in the 7- 
metric was studied by Herrera et al. 13]. The 7-metric 
is axially symmetric and static. Thus as in the case of 
JNW metric the existence of the Killing vectors & = dt 
and = Off, implies that the following two quantities are 
constants of motion 



E = -g a ptfU? = FU l 



r 2 - 2Mr 
F 



u<t> 



(31) 
(32) 



where U a = (U t ,U r ,U e 1 U' t >). We assume that the 
geodesies are confined to the equatorial plane and thus 
impose 9 = \. 

Using the condition for the normalization of velocity 
U a U a — —1 and equations ([3"Tjl and (f3"2"]) the equation 
for the radial motion can be written as 



Gr 2 + V 2 = E 2 



(33) 



where V can be thought of as an effective potential and 
it is given by the following expression 



L 2 F 2 



V 1 = F 



2Mr 



(34) 
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B. Circular geodesies 

As described in the previous section the conditions for 
the circular geodesic are given by 



jjr = V r 



o, 



(35) 



where the dotted quantity means derivative with respect 
to the affine parameter, i.e. the proper time. These again 
are equivalent to 



V = E , and 



dV 
dr 



0, 



(36) 



meaning that the effective potential must be equal to the 
energy and it must admit an extremum. 

For the circular orbits to be stable against radial per- 
turbations the effective potential must admit a minimum, 
thus satisfying the condition 



d 2 V 
dr 2 



> . 



d*V 
dr 2 



(37) 

< 0, 
= 0. 

Now imposing condition (|36l) on the effective poten- 
tial (|34p we obtain the following expressions for angular 
momentum and energy 



For the unstable circular geodesies we have 
whereas for the marginally stable case we get 0^ 



L 2 = 



E 2 = 1 - 



f 2M y 



r 2 Mj 



- M(l + 2 7 ) ' 
M(l + 7) 



r ) r-M(l + 2 7 ) 



(38) 
(39) 
2 7 ) is 



Circular orbits exist if the condition r > M(l 
satisfied which ensures that E and L are real. 

Conserved energy and angular momentum blow up at 
the radius of the photon sphere which is given by the 
equation 



r ph = M{l + 2 1 ) 



The condition for the stability of the circular geodesies 
from equations is then 



where 



r > r + or r < r_ 



r± = M (l + 37 ± V57 2 - 1 



(41) 



(42) 



are the radii of marginally stable orbits. 

In Fig. [7] are shown the values of radial coordinate for 
the singularity, photon sphere, r + and r_ as the param- 
eter 7 is varied. Three distinct ranges of the parame- 
ters emerge where the structure of the circular geodesies 
will be qualitatively different, namely 7 € (0, 7 € 

(-4=, |) and 7 € (|,oo). We discuss these cases sepa- 
rately. 



2.5 - 




FIG. 7: The radial curves (in units of ^7) for the singularity 
(x-axis), photon sphere (thick line), r+ (dotted and dashed 
line) and r_ (dashed line) respectively as functions of the 
parameter 7. It is obvious from the graphs that the structure 
of the circular geodesies will be qualitatively different in the 
three different regimes of the parameter, namely (0,-jg), (775; 
|) and (5,00). 
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(40) FIG. 8: The conserved energy per unit mass E is plotted 



against the radius of the circular orbit r (in units of 2i0")^ or 
different values of 7. The curves located higher in the diagram 
correspond to higher values of 7. For values of 76 (i, 00), .E 
diverges at the singularity, decreases to a minimum at a finite 
radius and then increases to diverge again as r goes to infinity. 
For values of 7 6 (t^j \)i ^ goes to zero at the singularity 
and it has a maximum and a minimum. Whereas for values 
of 7 £ (0, -jg), E is monotonically increasing from zero at the 
singularity to infinite as r goes to infinity. 



The conserved energy and angular momentum as func- 
tions of the radius of the orbit of the test particle are 
plotted in Figs. [8] and |9] respectively for different values 
of the parameter 7. The qualitative behaviour of these 
functions is also manifestly different in the three regimes 
of the parameter mentioned above. 



8 



L 



6 










5 










1 










3 
2 










1 





1 2 3 4 5 6 2M 

FIG. 9: The conserved angular momentum per unit mass L 
is plotted against the radius of the circular orbit r (in units 
of Trjg) for different values of 7. The curves located higher 
in the diagram correspond to higher values of 7. For values 
of 7 £ (i 00) L diverges at the singularity and as r goes to 
infinity, thus having a minimum at a finite radius. For values 
of 76 ( -7= , i ) L has both a minimum and a maximum, goes to 
zero as r approaches the singularity and diverges as r goes to 
infinity. Whereas for values of 7 € (0, -j-) L is monotonically 
increasing from zero at the singularity to infinity as r grows 
indefinitely. 

1. Case 7 6 (0, ^ ) 

In this range of parameter r+ and r_ do not take real 
values and we have 

^ sing ^* , (43) 

meaning also that the photon sphere does not exist. This 
implies that the stable orbits can exist from infinity to 
the singularity: 

r G (r sing ,oo) . (44) 

In Fig JTU] the general behaviour of the effective poten- 
tial for this range of parameter is shown. The effective 
potential admits only one minimum. Since the potential 
does not admit maximum, unstable circular orbits do not 
exist. This is consistent with the existence of stable cir- 
cular orbits everywhere as discussed above. 

The qualitative behaviour of the functions E(r) and 
L(r), as shown in Figs. [8] and |9j is decreasing as the 
radial coordinate decreases from infinity and vanishes at 
the singularity. 

2. Case 7 G (-^, ±) 

In this range of parameters we get the following rela- 
tion 

r+ > r_ > r sing > r ph , (45) 



v 
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FIG. 10: The effective potential V is plotted against radial 
coordinate r in units of gjjy for 7 = 0.3 S (0, 4= ) for different 
values of L. The effective potential admits only a minimum 
implying that unstable orbits in this case do not exist. 
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FIG. 11: The effective potential for radial motion V is plotted 
as a function of radial coordinate r in units of ^7 for 7 = 
0.47 G ("7?i h) f° r different values of L. It can be seen that 
the V admits two minima and one maximum indicating the 
presence of two regions of stable circular orbits separated by 
a region of instability. 

which means that no photon sphere is present in the 
spacetime and there are two distinct regions where cir- 
cular geodesies can exist, namely 

r G (r+,oo) and r € (r sing ,rJ) (46) 

where the first region extends from infinity and the sec- 
ond region extends until to the singularity. 

In Fig. [TTJ we show the general behaviour of the ef- 
fective potential for this range of parameter values. The 
effective potential admits two minima. These correspond 
to the two stable circular orbits in the two regions we dis- 
cussed above. There is also one local maximum between 
the two minima which lies in the unstable region between 
the two stable regions. 

A particle reaching the innermost stable circular orbit 
of the outer disk would plunge in towards the singularity 
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FIG. 12: The effective potential for radial motion V is plot- 
ted as a function of radial coordinate r (in units of ^rj) for 
7 = 1.4 £ (|,oo) for different values of L. It can be seen 
that the V admits one minimum and one maximum, which 
corresponds to stable and unstable circular orbit respectively. 



this range. The effective potential admits one maximum 
and one minimum which correspond to the unstable and 
stable circular orbits respectively. The unstable circular 
orbits lie in the region between the photon sphere and 
r+- 

The qualitative behaviour of the energy and angular 
momentum functions E(r) and L(r) shows a minimum 
and it diverges at a certain radius above the singularity 
corresponding to the photon sphere (see Figs. [5] and [5]). 
This is consistent with the above discussion on the struc- 
ture of geodesies since as r decreases from infinity, the 
stable circular geodesies can exist up to the minimum, 
then unstable orbits exist in the region between the min- 
imum and the radial value where E{r) and L(r) blow up, 
which is a photon sphere. No circular orbits can exist 
below the photon sphere. 



IV. ACCRETION DISK PROPERTIES 



until it reaches r_ where stable orbits are allowed again. 
It would then shock and circularize its motion spiralling 
towards the singularity. 

The qualitative behaviour of the energy and angular 
momentum functions E(r) and L(r) is shown in Figs. [5] 
and [5J In this case E and L vanish at the singularity 
and both have a local maximum and a minimum indicat- 
ing again that there are two regions where stable circular 
geodesies can exist. The two regions would correspond 
to the places where functions decrease with decreasing 
r, namely from infinity to first minimum and from the 
maximum to the singularity. Again the unstable region 
extends between maximum and minimum, where E(r) 
and L(r) increase with decreasing r. Since these func- 
tions do not blow up at any radius we infer that the 
photon sphere is absent in this case. 

3. Case 7 G (|, oo) 

For this range of values of the parameter we have the 
following relations 

r+ > r ph > max{r_,r sm9 } , (47) 

as can be seen from Fig. [7J Note that r_ can be cither 
larger or smaller than r S i ng but in both cases lies below 
the photon sphere. The stable circular orbits can exist 
in the region 

r e (r+,oo) , (48) 

which is above the photon sphere. Unstable circular or- 
bits can exist between the photon sphere and r + . No 
circular orbits can exist below the photon sphere until 
the singularity. 

In Fig. [T2] we show the general behaviour of the effec- 
tive potential as a function of radius for the parameter in 



In the context of astrophysics an accretion disk is a 
nearly flat disk of heated gas whose particles slowly spiral 
onto a central accreting compact object. Observationally 
such accretion disks can be detected by their emitted ra- 
diation and the spectrum of this radiation can be a valu- 
able source of information about the properties of the ac- 
cretion disk itself and the central compact object. In the 
present framework we approximated the accretion disk 
with an infmitesimally thin disk of test particles moving 
on circular orbits. The gas is then represented by the test 
particles moving with angular velocity u> = ^ on Kep- 
lerian orbits around the central compact object. Then 
E(r) and L(r) would represent the energy and angular 
momentum of the test particles circling on the r-orbit of 
the disk. Then the last stable circular orbits represent the 
limits within which the disk can exist with the innermost 
stable circular orbit given by the radius r + representing 
the inner edge of the outer disk. A test particle reaching 
r + would then plunge in free fall until it reaches either the 
singularity or another radius where stable circular orbits 
are allowed. The general relativistic setting for accretion 
disks around a black hole was developed by Novikov and 
Thorne [HI and by Page and Thorne [l5[ for steady state 
accretion on the equatorial plane onto a rotating black 
hole and it can be extended to the scenarios with naked 
singularities analyzed here. 

Page and Thorne assumed the spacetime to be station- 
ary, axially symmetric, asymptotically flat and reflection 
symmetric in an equatorial plane. Here both metrics sat- 
isfy these requirements once we consider accretion disks 
on the equatorial plane (note that the JNW metric is 
spherically symmetric and thus such a choice of the plane 
is arbitrary while the 7 spacetime is axially symmetric 
and therefore the equatorial plane is fixed). Further- 
more the choice of the radial coordinate complies with 
the requirements imposed in 15] as can be seen from the 
asymptotic behaviour of both metrics and from the fact 
that r reduces exactly to the radial coordinate of the 
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FIG. 13: Radiant energy flux for unit mass accretion rate 
in natural units for the JNW spacetime with fj, < i/E. Note 
that the lowest curve represents the flux for the Schwarzschild 
black hole. 

Schwarzschild metric for 7 = 1 in the 7-metric and to 
a translation of the Schwarzschild radial coordinate for 
(j, = 1 in the JNW metric. 

From the knowledge of the structure equations describ- 
ing the gas in the disk it is possible to construct the quan- 
tities relevant for observations [15| . The radiant energy 
flux can be written as 

where Mq is the mass accretion rate and it is assumed to 
be constant for steady state accretion. Note that since 
both metrics considered here are continuously linked to 
the Schwarzschild solution one retrieves precisely the for- 
mula for the radiant energy flux in the Schwarzschild 
spacetime once the values fj, = 1 or 7 = 1 are chosen. 

As an example we studied here the qualitative behav- 
ior of the radiant energy flux in the JNW metric and in 
the 7-spacctime in natural units (therefore to obtain the 
effective flux one would have to reintroduce c and G in 
the equation), for different values of the parameters /j, 
and 7. 

The general feature that emerges for the JNW metric 
is that when there is an inner boundary for the disk one 
obtains a flux qualitatively similar but of greater magni- 
tude with respect to that of a black hole (see Fig. [T5|) . 

On the other hand when the accretion disk is allowed 
to extend up to the singularity the total flux is diverging 
in the limit of r going to the singularity (see Fig. [Til) . 
This divergence of flux in the limit of approach to the 
singularity is clearly unphysical. Nevertheless we note 
here that as the outgoing radiation becomes larger, the 
outward force acting on the infalling accreting matter 
will also become bigger. It can be shown that at a par- 
ticular radius the pressure exerted due to the outgoing 
radiation will balance out the inward gravitational pull 
of the gas in the disk. Therefore the accretion disk will 
not extend below this radius regardless of the presence 



f/M„ 




FIG. 14: Radiant energy flux for unit mass accretion rate in 
natural units for the JNW spacetime with fi > V5- Note that 
the lowest curve represents the case fi = \/5, with all values 
greater than y/E diverging as r approaches the singularity. 



of a singularity at the center. Thus although the circular 
geodesies as we have shown in this paper can extend all 
the way up to the singularity, it is reasonable to expect 
that physically the accretion disk will be terminated at 
a finite radius above the singularity. Nevertheless this is 
a complex issue that is not entirely well understood at 
present and requires more work and further discussion 
which is beyond the scope of this paper. 

Furthermore since the 7-metric is not spherically sym- 
metric it should be noted that the axial symmetry will 
have some influence on the structure of accretion disks 
close to the singularity. In fact it can be shown that close 
to the center strong forces in the z direction appear in 
the spacetime. These forces would disrupt the inflow of 
particles on the equatorial plane thus deviating them to- 
wards the poles, which could be a mechanism that may in 
principle describe the formation of high energy jets from 
the poles of such compact objects. 

In the case of the 7-metric a similar behaviour is found 
for the radiant energy flux to that of the JNW case. The 
flux presents a behaviour qualitatively similar to that of 
a black hole when there is an inner edge (see Fig. IT5j) and 
diverges in the cases where the accretion disk extends all 
the way to the singularity (see Fig. H~6| . 

Note that since the total mass as seen by an observer at 
infinity for the 7-metric is M7 it can happen for values of 
7 < 1 that the flux is smaller than that of a Schwarzschild 
black hole with mass M. Nevertheless we expect it to 
be greater than that of a Schwarzschild black hole with 
mass Mj, which would be the right spherical source to 
compare it with. 

To summarize we studied in detail the structure of the 
circular geodesies in JNW and 7 spacetimes in the equa- 
torial plane. JNW and 7-metric are obtained by variation 
of the Schwarzschild black hole by addition of the mass- 
less scalar field and by deformation of spherical symme- 
try to oblate/prolate spheroidal geometry respectively. 
Naively speaking the event horizon in the Schwarzschild 



11 



f/M„ 
5-10' 4 -i 




FIG. 15: Radiant energy flux for unit mass accretion rate in 
natural units for the 7-spacetime with 7 > 1/2. 




23456789 10 



FIG. 16: Radiant energy flux for unit mass accretion rate in 
natural units for the 7-spacetime with 7 < 1/2. Note that 
the flux diverges for values of 7 < 1/yE. 

black hole transforms into a naked singularity. These so- 
lutions are described by two parameters, one related to 
the mass as in the Schwarzschild case and another related 
to the scalar charge (/i) or the deformation (7) and con- 
tinuously linked to the Schwarzschild spacetime. In fact 
the Schwarzschild metric is obtained by setting fi = 1 
and 7 = 1. The structure of the circular geodesies un- 
dergoes a qualitative change as we vary the values of the 
extra parameter. In both spacetimes the whole range of 
the parameters can be divided into three regimes. 

- For fi £ (1, 2) in JNW metric and for 7 £ (|, 00) in 
7-metric a photon sphere surrounds the singularity 
and stable circular geodesies can exist from a ra- 
dius above the photon sphere to infinity. Unstable 
circular geodesies can exist below this radius up to 
the photon sphere. No circular geodesies can exist 
below the photon sphere. 

- For fj, £ (2, in JNW metric and for 7 £ (-4=, ±) 
in the 7-metric, stable circular geodesies can ex- 
ist from the singularity to a finite radius and from 



a larger radius to infinity. There is a region con- 
taining unstable circular geodesies between the two 
stable regions. No photon sphere is present. 

- For n £ ("v/5, 00) in JNW metric and for 7 £ (0, 
in 7-metric, the circular geodesies exist everywhere 
from the singularity to infinity. Again no photon 
sphere is present. 

This study suggests that there is a significant difference 
as far as the structure of circular geodesies is concerned 
in JNW and 7 spacetime containing naked singularity as 
compared to the Schwarzschild black hole, which would 
reflect into their accretion disk properties. 

Similar analysis of circular geodesies in connection 
with properties of accretion disks around naked singu- 
larities have recently drawn some attention as we try to 
understand the nature of compact objects such as the 
supermassive ones that exist at the center of galaxies. 
Harko and Kovacs Q studied a modified Kerr metric with 
scalar field where, due to the presence of the scalar field, 
a naked singularity is present in the spacetime. They 
showed that for certain ranges of the parameters, stable 
circular orbits can exist up to the singularity, a result 
that does not hold for the Kerr black hole (obtained in 
the limit of no field) . 

A similar analysis has been carried out for the Kerr and 
Reissner-Nordstrom metrics by Pugliese et al @ where it 
was shown that for the Kerr naked singularity (a/M > 1) 
unstable circular orbits can extend up to the singularity 
in some cases, though this does not happen for stable 
circular orbits. Therefore an accretion disk made of test 
particles always extends from a minimum radius to in- 
finity, without any forbidden regions or regions of insta- 
bility. The disk can be divided into substructures based 
on whether they consist of corotating and/or counter- 
rotating particles. For a naked singularity, depending 
on the value of a/M, the accretion disk may be divided 
into: (a) two parts - an inner disk with co-rotating parti- 
cles and an outer disk with both co-rotating and counter- 
rotating particles, (b) three parts, namely an innermost 
disk of counter-rotating particles, a second disk of co- 
rotating particles and an external disk with both. For a 
Kerr black hole the allowed region of circular orbits, sta- 
ble or unstable, terminates outside the horizon. However 
the accretion disk is similar, extending from a minimum 
radius to infinity and it can have only two such substruc- 
tures, a co-rotating inner disk and an external disk with 
both kinds of particles; here the size of the inner disk has 
an upper limit and can even go to zero. In the case (a) for 
naked singularity, the inner disk has a minimum size but 
there is no upper limit thus showing that there is some 
difference in the accretion disk structure between Kerr 
black hole and Kerr naked singularity. In the Reissner- 
Nordstrom spacetime a naked singularity is obtained for 
Q/M > 1, where Q is the charge parameter. In this 
case there exists a forbidden zone outside the horizon or 
the singularity where no circular orbits are allowed. For 
the black hole case, the stability region and hence the 



12 



accretion disk is continuous from a minimum radius to 
infinity. In the case of a naked singularity, there can be 
two scenarios depending on the value of Q/M, namely 
(a) an accretion disk extending from a minimum radius 
to infinity, or (b) an accretion disk comprising an inner 
disk and an outer disk separated by a zone of instability. 
Therefore the accretion disk structure around a Reissner- 
Nordstrom naked singularity shows some difference with 
respect to the black hole case only for a certain range of 
the parameter. 

V. CONCLUDING REMARKS 

The study of circular geodesies provides only the first 
tool to model accretion disks around a compact source. 
In order to obtain measurable quantities that could even- 
tually be checked against observations, a detailed analy- 
sis of the efficiency of conversion of mass of the infalling 
particles into emitted radiation and the luminosity spec- 
trum of radiation of the disk are needed. For example 
assuming that the accretion disk is in thermodynamical 
equilibrium one can approximate the radiation emitted 
by the disk with a black body spectrum where the tem- 
perature T of the black body is related to the energy flux 
by f{r) = crT 4 (r) (with a being the Stefan-Boltzmann 



constant). Then, considering redshift and inclination of 
the disk with respect to the observer, it is possible to 
evaluate the observed luminosity spectrum of the black 
body, which is an observable quantity Q- 

It is reasonable to suppose that the pressure of the out- 
going radiation near the singularity will be very strong, 
thus balancing the inflow of particles from the accretion 
disk. This effect could possibly drive the particles away 
from the equatorial plane and it should also be analyzed 
in some detail. Furthermore electromagnetic fields and 
possible perpendicular forces that could deviate the par- 
ticles towards the poles thus creating high energy jets 
could also be considered. 

Nonetheless the analysis such as given here makes it 
clear that if naked singularities which are hypotheti- 
cal astrophysical objects do exist in the universe, they 
would bear an observational signature considerably dif- 
ferent from that of a black hole of the same mass. There- 
fore it is quite possible that in the foreseeable future we 
will be able to devise some observational tests to check if 
such exotic objects are indeed present somewhere in the 
universe. 
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